Dy2Ti207 Spin Ice: a Test Case for Emergent Clusters in a Frustrated Magnet 
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Dy2Ti207 is a geometrically frustrated magnetic material with a strongly correlated spin ice 
regime that extends from 1 K down to as low as 60 mK. The diffuse elastic neutron scattering 
intensities in the spin ice regime can be remarkably well described by a phenomenological model 
of weakly interacting hexagonal spin clusters, as invoked in other geometrically frustrated magnets. 
We present a highly refined microscopic theory of Dy2Ti207 that includes long range dipolar and 
exchange interactions to third nearest neighbors and which demonstrates that the clusters are purely 
fictitious in this material. The seeming emergence of composite spin clusters and their associated 
scattering pattern is instead an indicator of fine-tuning of ancillary correlations within a strongly 
correlated state. 



Geometrically frustrated magnetic materials have in 
recent years furnished many new paradigms for the explo- 
ration of novel condensed states of matter. Examples in- 
clude spin ice behavior [J, [2, i, i, i, 0, 0, i, i, , col- 
lective paramagnetism or spin liquid behaviour |l2l| , spin 
glassiness in stoichiometrically pure materials [13| and an 
anomalous Hall effect . Frustration is detrimental to 
the development of conventional periodic arrangements 
of magnetic moments at low temperatures [l5l |. As a 
consequence, at temperatures that are small compared 
to the scale of the leading interactions, highly frustrated 
systems typically form collective dynamical states [l6| . 
Most commonly, the massive degeneracy of such "cooper- 
ative paramagnetic" states is ultimately resolved by spe- 
cific material-dependent perturbative effects. This is the 
origin of the diversity of experimentally observed phe- 
nomena [isj . 

Among the many systems that exhibit geometric frus- 
tration [l5| , all the above phenomena have been observed 
in materials with magnetic moments (spins) residing on 
a network of corner sharing tetrahedra, or pyrochlorc lat- 
tice. This arrangement is found, for example, in the rare 
earth titanatcs i?2Ti207; with magnetic rare earth ions 
(i?^+ = Ho, Dy . . . ), and in the spinels hke AB2X4 {A = 
Zn, Cd) , with magnetic transition metal ions (B = Fc^"*" , 
Cr'^+). In this paper we are concerned with Dy2Ti207, 
a spin ice [H, 0, ^ . In this material, the local spin corre- 
lations arc characterized by the ice rule: two spins point 
in and two spins point out of every tetrahedron [J, |3, ■ 
The strongly correlated spin ice regime (analogous to the 
spin liquid or collective paramagnetic regime in geometri- 
cally frustrated antiferromagnets [31) extends from 1 K 
down to the lowest measured temperatures (< 100 mK). 

Given the many phenomenologies exhibited by frus- 
trated materials, the recent idea that there may 
be an organizing principle that describes the spin cor- 
relations in spin liquid states is very appealing. Us- 
ing inelastic neutron scattering data, it was shown in 



Ref. [17[ that the spatial correlations of magnetic exci- 
tations in the spinel ZnCr2 04 are well described by a 
model of strongly bound hexagonal spin clusters that 
are uncorrelated with respect to each other. Subse- 
quently, similar evidence of such clusters was obtained 
in CdFe204 [ll, CdCr204 [H and Yo.5Cao.5BaCo407 
[20! . lending weight to the idea that the clusters may ef- 
fectively be emergent objects. In most cases it is numer- 
ically or analytically rather intractable to examine these 
systems by detailed microscopic calculations. However, 
Dy2Ti2 07 [31 provides an exception, as we show here 21 1. 

In Ref. [7|, we and collaborators showed that the 
well-established microscopic "standard dmolar spin ice 
model" (s-DSM, defined below 0, ^ U^) describes 
the spin correlations in the spin ice state of Dy2Ti207 
with only modest success. The original objective for 
the present study was to understand the reason for the 
differences between theory and experiment. To do so, 
we examined two models of the spin-spin correlations in 
Dy2Ti207 by calculating the Fourier image, /(q), of the 
spin correlation function and comparing it to the energy- 
integrated neutron scattering structure factor 0]. The 
Dy^+ ion in Dy2Ti207 has a Kramers doublet as single- 
ion crystal field ground state [2J| that, for each site i, is 
well approximated by a classical Ising degree of freedom, 
Si = ±1, defined along its local [111] trigonal axis Zi [i^ ]. 
We calculate /(q) using 
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where (. . .) denote thermally averaged (siSj) correlations 
between the Ising spins at sites i, j; z^ is the component 
of the quantization direction at site i perpendicular to the 
scattering vector q, N is the number of spins and /(|q|) 
is the Dy'^^ magnetic form factor [2^. For comparison 
with experiment, /(q) is adjusted by an overall scale fac- 
tor and a slowly varying linear-in-|q| background. The 
experimental diffuse neutron scattering, measured in the 
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elastic approximation in the (hhl) plane (from Ref. is 
shown in Fig.[lji. As previously observed Q, in addition 
to structural Bragg peaks (e.g. (004), (222)) and bright 
broad features at (001), (003) and (f ||), the experimen- 
tal /(q) is further decorated by hexagonal loops of diffuse 
scattering running along the Brillouin zone boundaries. 

The first model of correlations we examine is phe- 
nomenological and based on postulating an ansatz for 
{siSj) in (1). We assume that (siSj) can be viewed as 
generated by static clusters that remain uncorrelated be- 
tween themselves. Each cluster is a zero-magnetization 
hexagonal loop of spins that circulate perpendicular to 
the loop normal. These clusters are the discrete equiva- 
lent of the "emergent" clusters used to describe the in- 
elastic /(q) in ZnCr204 [H (see Fig. Ic in Ref. 
Taking into account that all spins on a pyrochlore lattice 
can be grouped into non-overlapping hexagons without 
breaking the ice rules, that hexagon normals have four 
possible orientations, and that each hexagon has two pos- 
sible senses of "spin circulation" around the normal, /(q) 
can then be calculated using Eq. (1). Fig.[T)D shows /(q) 
calculated using the spin cluster scattering function. This 
model describes the experimental data quite well. The 
selection of hexagonal clusters as effective degrees of free- 
dom in Dy2Ti207 does not incorporate any microscopic 
information about the host material, and thus be viewed 



as another example [17| of emergent composite spin clus- 
ters in frustrated systems. 

The second approach we use to determine (siSj) is mi- 
croscopic. The s-DSM was previously shown to account 
fairly well for the spin ice phenomenology of Dy2Ti2 07 
0, [H, [2^. It comprises the magneto-static dipole in- 
teraction, which gives a ferromagnetic nearest neighbor 
coupling, that competes with a weaker antiferromagnetic 
nearest- neighbor exchange interaction Ji. Sufficiently 
strong antiferromagnetic Ji would lead to long range 



order [22, |23[. It is the interplay between the prop- 



erties of the spin ice manifold and the symmetry and 
long range nature of the dipolar interaction that leads 
to a correlated spin ice state over an extended tempera- 
ture range 22, 23, 27 1. Surprisingly, the s-DSM is much 
less successful than the simple phenomenological cluster 
model at describing /(q). 

Monte Carlo (MC) simulations have shown [tI that the 
s-DSM correctly describes the location and relative inten- 
sity of the strong /(q) features, but fails to reproduce the 
hexagonal zone boundary scattering (ZBS), Fig. [TJ;. We 
interpret this as a sign that the s-DSM is incomplete and 
needs to be extended [111, Ell ■ 

Dy2Ti207 displays a number of phase transitions and 
other structured resp onse in an applied magnetic field 
H [3, i, i, S, S, i, [Ifl El • Independently of the disagree- 
ment between Fig.[T^ and Fig. [it, the necessity to adjust 
the s-DSM was previously also suggested by the obser- 
vation that, while it qualitatively explains the in-field 
transitions, the s-DSM does not correctly predict their 
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FIG. 1: (color online). Neutron scattering intensity /(q) of 
Dy2Ti207 as a function of the wavevector. (a) Experimental 
elastic scattering intensity at 300 mK 0] in the (hhl) planes 
of the reciprocal space is notably amassed along the hexago- 
nal zone boundaries, (b) /(q) calculated using the model of 
uncorrelated hexagonal spin clusters (see text), (c) /(q) ob- 
tained from MC simulations on the s-DSM of Dy2Ti207 
describes the main experimental features, but is inadequate 
in reproducing the ZBS (d) The g-DSM allows for an ex- 
cellent match between the theoretical and experimental /(q), 
and allows us to identify that correlations beyond 3rd nearest 
neighbors are the microscopic origin behind the ZBS. (e,f) 
Quantitative comparison of the experimental (a) and theo- 
retical (b, c, d) data along a cut through the reciprocal space 
chosen to emphasize the mismatch with the s-DSM. Panels c, 
d, e, f produced by MC simulations of 8192 = 8"^ x 16 spins. 



temperatures [28|, unless properly adjusted by pcrturba- 
tive exchange couplings beyond Ji. Here, we consider a 
generalized dipolar spin ice model (g-DSM) of Dy2Ti207 
that includes second J2 and third J3 nearest neighbor 
exchange couplings: 



i>j 
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Here, i, j span the sites of the Dy'^+ ions, rtj and fij 
is the length and direction of the vector separation be- 



3 




-0.4 

3.1 3.2 3.3 3.4 3.5 
J,,K 



0.41 

3.1 3.2 3.3 3.4 3.5 



0.20 



0.15 



0.10 



10.05 
0.00 



FIG. 2: (color online), rms deviation a of Monte Carlo 
vs. experimental Q specific heat data in the Ji — J2 plane 
at fixed value of J3 = 0.025 K. The isolines are drawn at 
o- = 0.05,0.10,0.15 and 0.20 J moP^K"^ outward. The two 
functions a [(a), obtained for the ri=0.7 K-1.49 K tempera- 
ture interval) and (b), for r2=1.5 K-2.8 K], yield candidate 
solutions (delineated by inner isolines) in overlapping parts of 
the Ji — J2 plane. The determination of the couplings in 
the Hamiltonian ([2]) of Dy2Ti207 is obtained via a fitting to 
several bulk experimental thermodynamic data, akin to the 
procedure illustrated here (see text). 



tween spin pairs. D = ^^{{J'^gijL^Y / {A-Kri^) = 1.3224 
K is the strength of the dipolar interaction at nearest- 
neighbor distance n = a'\/2/4, as obtained from the esti- 
mate {J') = 7.40 [29| within the ground state doublet of 
Dy3+ in DyzTiaOr, and from the size a = 10.124 A ^ 
of the cubic unit cell of the material. 5 = 4/3 is the T)y^'^ 
Landc factor, is the Bohr magneton, and is the ex- 
change coupling of spins at distance r^. In reality, there 
should be two different 3rd nearest neighbor exchange 
couplings. We take them equal, which, as shown below, 
provides a consistent description of Dy2Ti207. By do- 
ing mean-field theory (MFT) calculations of J(q) in the 
paramagnetic regime [30[ we checked that the conclusions 
about /(q) are not sensitive to ^ 50% deviations of the 
two J3 from equality. 

Rather extensive MC simulations were used to investi- 
gate whether the g-DSM can provide a consistent descrip- 
tion of the phenomenology of D^Ti207 both for H = 
[3, i, d and H 7^ i6i, J, i, [1^ regimes. MC simula- 
tions were performed on system sizes of 1024 = 4^ x 16 
spins with periodic boundary conditions, treating the 
long range dipolar interaction by the Ewald method [2^ . 
For each set of parameters in ^ and each temperature, 
10^ MC steps per spin for equilibration followed by an 
additional 10^ steps per spin for production were per- 
formed. The simulations exploited both single spin- flip 
and loop-update [2^ Metropolis algorithms, the latter 
being a generalized algorithmic MC version of the dy- 
namics of the hexagonal cluster degrees of freedom dis- 
cussed above (cf 



for H nearly along the [112] direction [ol. [28| . F or fixed 
_D, the transition is controlled by J3 only [28|. Map- 
ping the MC critical temperature, obtained as location 
of the maximum in magnetic heat capacity Cm, to the 
experimental estimates of 0.34, 0.28(1) and 0.26(1) K, 
obtained b y sp ecific heat 0, susceptibility [§| and mag- 
netization [10| measurements, we conclude that 0.019 K 
^ ^3 £ 0.026 K. Second, we examine Cm{T) data in 
H = 0, 0, Si for constraining Ji, J2, as C,„(T) is 
only weakly sensitive to the allowed small variation of 
J3. We determine optimal Ji, J2 by minimizing the root 
mean square (RMS) deviation a of experimental vs the- 
oretical temperature curves for Cm(T)/T. The function 
cr( Ji, J2) at J3 = 0.025 K is plotted in Fig.[2]using Ref. Q 
data. The two panels illustrate that a{Ji , J2) has a stable 
minimum even if determined over two different temper- 
ature intervals. Third, we exploit an empirical equation 
He = 0.90(1) + O.OST, T < 0.36 K @, i] for the line of 
phase transitions in Dy2Ti207 for H along [111]. We ob- 
tain values of He at several temperatures as positions of 
maxima on MC C„i(H) curves, and vary Ji and J2 to 
match the experimental Hc(r). This, the above H = 
analysis and an additional observation that at H = 
a magnetic phase transition, if any, can occur only at 
Te < 300 mK 0, 0, cf. Fig. [TJi, yields: 3.26 K < Ji < 
3.53 K; —0.20 K < J2 < K. Our value for J2 is consistent 



with the value J2 = —0.1 K reported in Ref. [ll| from 
fitting© with D = 1.41 K, Ji = 3.72 K, J3 = 0.03 K 
(Ref. [28[) to the experimental Cm(H) for H along [111]. 

To verify whether the experimental /(q) in Fig. [T^ 
can be described by we perform MC simulations 
of the optimized g-DSM at 300 mK. We find that the 
g-DSM with Ji, couplings within the above allowed lim- 
its is consistent with the experimental scattering pattern 
in Fig. [1^, but requires even stronger restriction on J2: 
-0.16 K < J2 < -0.10 K. For instance, the theoreti- 
cal pattern at Ji = 3.41 K, J2 = -0.14K, J3 = 0.025K 
(Figs. [TJi, f) reproduces the experimental one (Fig. [1^) 
extremely well. As seen from Figs. [TJ;, d, e and f, 
the small adjustments to the initial s-DSM [S^] do re- 
sult in a redistribution of the scattering response, with 
intensities that now capture correctly both major fea- 
tures of the experimental pattern, and the weaker ZBS. 
The characteristic features on the experimental /(q) for 
Dy2Ti207, first described above as arising from hexag- 
onal spin clusters, can therefore be "straightforwardly" 
described within a conventional microscopic treatment of 
the static spin-spin correlation function [31 ] . 

Having now a credible microscopic model in hand, are 
we able to explain the success of the phenomenological fit 
(Fig. [TJj)? To answer this question, we examine the di- 
Fig. 1 [23|). Finite size effects were rect space {siSj) correlations behind the reciprocal space 

patterns in Fig. [Th), c, d. We intimate that the picture of 
independent clusters is equivalent to a correlation func- 
tion truncated beyond third nearest neighbor distance 
(i.e. outside the hexagonal cluster). The (siSj) obtained 



found to be insignificant for the discussion below. 

Three groups of experimental data about Dy2Ti2 07 
are used to determine the in ([2]). First, we consider the 
temperature-driven ferromagnetic ordering of Dy2Ti2 07 
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from MC simulations on both the s-DSM (Fig. [T]:) and 
optimized g-DSM (Fig. [TJi) are very similar at short dis- 
tances. If truncated beyond third nearest neighbor dis- 
tance, the Fourier transform of each produces a pattern 
close to that of experiment (Fig. [T^) . We conclude that 
the difference between the two MC /(q) patterns (Fig.[TJ; 
and Fig. [TJi) is caused by correlations beyond third neigh- 
bors. For the tuned g-DSM, the long range part strongly 
reinforces the ZBS arising from the short range part, 
contributing the majority of this scattering. This con- 
tradicts the notion of independent hexagons, which un- 
derly the calculation of /(q) in Fig.[T|3, since correlations 
between hexagons are of the same order as the correla- 
tions defining the hexagons. Thus, weakly interacting 
clusters are not an appropriate microscopic description 
of Dy2Ti207, nor are they here an effective organizing 
principle. Rather, the ZBS are caused by subsidiary 
{ J2 and J3) interactions that are largely inconsequen- 
tial to the formation of the strongly correlated spin ice 
regime 22, l23|, |27| . We further confirmed the fine-tuned 
nature of the ZBS using MFT to calculate /(q) in the 
paramagnetic regime [30| . Finally, we note that indepen- 
dent hexagonal clusters, unlike the s-DSM and g-DSM, 
fail to predict the "pinch point" scattering characteristic 
of a spin ice manifold, which is disccrnablc in the exper- 
iment. 

We finally comment on the broader consequences that 
the observed complementarity of the phenomenological 
and microscopic interpretations of neutron scattering 
patterns might have for other systems. Our results sup- 
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port the idea expressed in Ref. [l7| that cluster-like scat- 
tering (CLS) may be common among highly frustrated 
magnets and establish here a microscopic mechanism by 
which it can arise. However, by considering Dy2Ti207 as 
a test case, we have shown that the CLS does not neces- 
sarily imply the emergence of "real" clusters [s^l, nor a 
new organizing principle. Instead, the CLS is the prop- 
erty of a strongly correlated liquid-like state and a con- 
sequence of the sensitivity of frustrated systems to per- 
turbations [Til . Our study also shows how neutron scat- 
tering can allow one to fully interpret CLS provided suf- 
ficiently accurate data and a computationally tractable 
theoretical microscopic model are available. 
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